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Abstract 

We argue that the geodesic rule, for global defects, is a consequence of the randomness 
of the values of the Goldstone field in each causally connected volume. As these vol- 
umes collide and coalescence, (p evolves by performing a random walk on the vacuum 
manifold Ai. We derive a Fokker-Planck equation that describes the continuum limit of 
this process. Its fundamental solution is the heat kernel on A4, whose leading asymptotic 
behavior establishes the geodesic rule. 
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I. INTRODUCTION 



To predict the density of topological defects formed after a phase transition, Kibble [1] 
postulated two "rules" : the randomness and independence of the value of the Goldstone 
scalars 0, parametrize the vacuum manifold Ai within each causally connected region, 
and the "geodesic rule" . According to the latter, when two regions of Ai having different 
values of (p collide with each other, interpolates between these values across the common 
boundary, following the shortest path [1]. The geodesic rule is simple, and as such, aes- 
thetically appealing and favorable for numerical implementation, but it lacks a thorough 
justification (see however [2], [3]). Moreover, there is no consensus on its validity, or even 
on its extensions in the case of local defect formation [4]-[8]. 

We argue that the geodesic rule is a result of the stochasticity of the selection of 
within each causal horizon. Our argument relies on the random nature of the collisions 
between causally unrelated regions on A4. It differs from the other treatments on this 
issue, which rely on energetic considerations [2], [3]. The present discussion applies, as is, 
to theories with broken global symmetries only. The presence of gauge symmetries would 
introduce several subtleties in the present arguments. These diflicultics are directly re- 
lated to the fundamental lack of understanding of the vacuum structure of gauge theories. 
Explicit understanding of the topology and geometry as well as paramctrizations of the 
underlying moduli space of the affine connections of gauge theories is still lacking in most 
cases of physical interest. As a result, even the formulation of appropriate stochastic 
equations, analogous to the Fokker-Planck, becomes problematic. For these reasons we 
will deal only with the case of global defects in this treatment. 

The organization of this paper is as follows: Section II presents the derivation of the 
geodesic rule and the assumptions on which this derivation relies. Section III analyzes 
the domain of applicability of this result and re-examines its assumptions. Section IV 
presents a brief discussion and our conclusions. 

II. DERIVATION OF THE GEODESIC RULE 
At the outset of our treatment, it is worth pointing out that the Goldstone field (order 
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parameter) is scalar only as far as its space-time properties are concerned. In a general 
classical field theory, is a section of a vector bundle V over a space-time M which is most 
frequently considered to be M'*. The typical fiber of P is a hnear representation, in most 
cases, of a Lie group G [9] . Then provides a local parametrization of G in V. In a quan- 
tum or thermal field theory we start with an appropriate space of maps / : M — > G, or, in 
this language, with an appropriate space of sections of V. These sections must have some 
physically acceptable regularity properties, most frequently pertaining to their differen- 
tiability, as well as appropriate normalization properties with respect to the symmetric 
or hermitian inner product of G and the Riemannian inner product of M [9]. Moreover, 
in order to ensure the finiteness of the energy of the system, we impose some asymptotic 
requirements on / which amount, very frequently, to effectively considering a compactifi- 
cation of M. Such maps tend to form, very frequently, Sobolev spaces on the sections of V. 

A classical/non-thermal field theory can be determined by its Lagrangian C [9], [10]. A 
quantum/thermal field theory can be defined by a path integral on appropriately defined 
functional spaces, as was mentioned above. The result of the calculation of such a path 
integral can be encoded, under some rather general conditions, in such a way as to give 
rise to extra terms to the classical potential Vc- The requirement of renormalizability 
[9], [10], when applicable, often imposes additional constraints on the nature of these ex- 
tra terms, thus severely constraining the form of the initially allowed form of C Our 
arguments rely on the existence of such an effective potential V^ff [10] for the description 
of the thermal field theory. Renormalizability constraints of the underlying theory are of 
no concern here, since we are dealing with effective, as opposed to fundamental, theories 
from the outset. Then, the effective potential allows us to describe the system in an ap- 
parently non-thermal manner and that is why many of our arguments have a very strong 
classical/non-thermal flavor. 

For concreteness, we begin by considering a first order phase transition proceeding by 
bubble nucleation [11]. This hne of argument essentially holds also for continuous tran- 
sitions proceeding by spinodal decomposition, as we comment at the end of the paper. 
We assume that the phase transition giving rise to the topological defects proceeds by 
thermalization and will be eventually completed. It is known [12], [13] that after their 
formation, bubbles whose radius is larger than a critical value Tcrit tend to expand and 
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bubbles whose radius is smaller than rcrit tend to shrink until they disappear. For this 
reason, we restrict our attention to bubbles of radii greater than rcrit- We are interested 
in length scales much bigger than rcrit- Then without loss of generality we can pretend 
that r can take any positive value with lower bound zero. 

A first assumption is that when two bubbles collide and coalesce, the bubble which is 
generated has an effective value of which is the weighted average of the values of the 
two bubbles that coUided. In this weighted average, the contribution of each of the two 
colliding bubbles should be proportional to the number of points it contains, i.e. to its 
volume. It is simplicity that makes such a choice plausible. It is very possible that there 
are systems for which such a rule for calculating the weighted average does not apply, 
although we fail to see a generic reason why it should not. Thus, if the two coalescing 
bubbles have values of Goldstone fields 0i and (j)2 respectively, the resulting bubble will 
have an effective value of given by 

1401 + 1^ (1) 

where Vi and V2 are the corresponding volumes of each bubble. 

A second assumption is that at least two distinct, widely separated time scales exist 
for such a system: tc which describes the average time between two consecutive bubble 
collisions and tm which is the mixing time of the two values of the Goldstone fields, after 
the bubble collision has taken place. These time scales are determined by the dynamics 
of each model. Their specific values, important as they are in the dynamics of phase 
transitions, are not relevant in our argument, so we will not delve into the issue of their 
calculation any further. This assumption is very strong as we expect that most systems 
violate it, more or less. However, as we explain in section III, for systems that follow 
it or weakly violate it, the geodesic rule should still hold, due to the wide domain of 
applicability of the Central Limit theorem. If tq -C tm, the values of (1) of the two bubbles 
will not have enough time to mix well between two collisions. This is a case favorable for 
the formation of topological defects [14], [15]. The opposite case, namely tq ^ tm, is very 
unfavorable for the formation of defects. For this reason we examine systems that belong 
to the former case. Incidentally, the probability of simultaneous collisions of three or more 
bubbles is so small that the contribution of such events to the number of produced defects 
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is negligible. 

Due to the reasons given in the first two paragraphs of this section, we need to describe 
aspects of the dynamics of the phase transition as it evolves through bubble collisions on 
Al X R. This evolution is parametrized by a scalar variable t e R. The relation between t 
and the actual time in which the phase transition takes place depends on the parameters 
of the specific model. We do not need to know any further details about such a relation 
for our argument's sake; it just suffices to know that such a relation exists. Because of 
the random distribution of bubbles and the random way in which they collide, (f) under- 
goes a random evolution i.e. it undergoes a random walk on M.. Let P{(l),t) denote the 
probability density function describing this stochastic process. To determine P{(l),t), we 
start by discretizing R, i.e. reducing it to a Euclidean one dimensional lattice, isomorphic 
to Z. We set the spacing of this lattice to be equal to tc since this is the shortest scale of 
t associated with the evolution of bubbles. Since tm ^ tc, we can assume without loss 
of generahty that tm is an integral multiple of tc- Then P{(t),t) has domain At x Z. To 
explicitly indicate this discretization, we rewrite t is as where i„ e Z. 

A third assumption in our argument is that the evolution of the system of bubbles 
described by P{(j), tn) is a Markovian process. This is actually an oversimplification. We 
have already assumed that tm tc, so after each collision the values of do not mix very 
strongly in the resulting bubble. Then the relative orientation of the bubbles before the 
collision took place leaves its trace on the resulting bubble through the distribution of the 
values of cj) on it. However, even within our approximation, the relative orientation of two 
bubbles before each collision is random. Eventually, we average over all possible coUisions 
between pairs of bubbles. The averaging process erases all trace of the relative orientation 
of bubbles before a coUision. In short, although the collision between two specific bubbles 
would give rise to a non-Markovian process, the randomness of orientation and the aver- 
aging process reduce the stochastic process to a Markovian one. This stochastic process 
is determined by P(0, i„) which obeys the Chapman- Kolmogorov equation [16], [20] 

Pih,^)^ [ P{4>^M4>2M) P{4>2M(t>iM) (2) 
Jm 

where P{(f)2, ^2|0i, ^i) denotes the conditional probability density for the stochastic variable 
to have value 02 at t2 given that it had value 0i at ti. The definition of the probability 
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density function P(0, tn) also gives 

P{4>, tn + Tc)= [ P{(l), tn + rcl^, t^) P{lp, t^) C^^ (3) 
JM 

We are interested in the behavior of the system at large t scales when compared to tc, as 
the system approaches its equilibrium, so we have to consider the continuum time limit 
Tc — > 0. Then t„ is replaced by t and a Taylor series expansion with respect to tc 
gives the master equation [16] 

= {wimni^. t) - w{mpi<t>, m # (4) 

Here denotes the transition probabihty per unit t for the Goldstone field to 

attain the final value ip, given its initial value was (f). 

A considerable simplification occurs, if we further assume that each bubble collides 
only with much smaller ones. This is a very strong assumption, which we expect to be vi- 
olated in may cases. Clearly the way the bubbles evolve and coUide depends on the details 
of the nucleation model describing the bubbles. It is moreover true that generally many 
bubbles of similar sizes collide with each other, a fact which renders this approximation 
invalid. The assumption is necessary though, in order to make the quantitative analysis 
tractable. We expect that because of the very wide range of applicabihty of the Central 
Limit theorem, the final conclusion that we reach applies to many cases that initially may 
even violate this assumption. We perform a Taylor series expansion in terms of the short 
displacement s that the scalar field experiences after each collision. Complications 
arise because is a Riemannian manifold which is not flat. The equations come as close 
as possible to the ones of the flat case by choosing a normal coordinate system. In this co- 
ordinate system the connection coefficients (Christoffel symbols) are zero at the origin O 
and, to a first order approximation the metric there is Euclidean. Let {ci, 62, ... , Cdim ju} 
be a coordinate basis of T^Ad. The master equation (4), when expanded up to quadratic 
terms in s in this basis [20], [21] gives 

= -V, [a\cP)P{4>,t)] + lv,V, [a^^{cP)P{ct>,t)] (5) 

where 

+00 

o^J(0) = J s's^ W{(t>;s) ds (6) 
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with 

W{i;\<f>)^W{<l>;s) (7) 

where s is the distance between and ip with respect to the metric of A4, and V is 
a Riemannian connection on M., not necessarily compatible with the metric [19]. This 
Fokker-Planck equation (5) has the first moment a'-{(j)) equal to zero since there are no 
"external fields", i.e. preferential directions on M.. Even if that is not the case, and the 
hnear term is non-zero, the stochasticity of the directions of the bubble collisions would 
eliminate such a linear term upon averaging over ToA4. Then (5) reduces to the diffusion 
equation 

dt 

with D'^{(f)) = i a^^(0) 



^V,Vj[D^\(l>)P{cP,t)] (8) 



If a bubble under consideration collides with either similar in size or much bigger 
bubbles, then the approximation leading to (5) is insufficient to describe the phenomenon. 
In such a case we need to keep more terms in the Kramers- Moyal expansion 

E „.,Jr'';,.. v..v.....v,...K.-»--->-x0)m^ (9) 

of the master equation (4). Actually, according to Pawula's theorem [20], which holds for 
a stochastic process in Euclidean space, the positivity of the transition probability rate, 
W{'ip\(f)), implies that the Kramers-Moyal expansion either terminates at the second order 
term (resulting in the Fokker-Planck equation (5)) or never terminates. Considering this 
theorem, our approximation leading to the termination of the Kramers-Moyal expansion 
(9) at the second order term is also optimal, in a sense, for technical reasons: any other 
choice would result in infinite sub-harmonic terms which would not be manageable without 
further, model-dependent, approximations. It is worth mentioning at this point that the 
expansion (9) is used heuristically. The covariant derivatives on A4 do not commute [19], 
since 

ViVj ek - Vj Vi Cfe = R{ei, ej) (10) 

where R{ei,ej)ek are the Riemann tensor components in the basis of ToA4 that we are 
using. The ambiguity in the order of the covariant derivatives appearing in the terms of 
order higher than two in (9), result in extra terms containing the Riemann tensor and 
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its covariant derivatives. The second order term may, at worst, acquire an extra term 
involving the scalar curvature of Al, a fact that does not alter the leading dependence of 
the resulting solution (13) on the distance function d{(f)), as will be seen in the sequel. The 
order, as well as the number of times, that the "diffusion" tensor components a*i*2 ' *'="'((/>) 
have to be differentiated depend on the details of the system. These ambiguities should 
be considered as an addition to the, essentially arbitrary at the mesoscopic level, choice of 
the Ito or the Stratonovich way of performing the stochastic integration [16], [21], which 
arises in the integration of the Fokker-Planck equation (5). In either case, a more detailed 
knowledge of the microscopic behavior of the system is required [16], [20] for the correct 
choices to be made. Such a particular choice, however, plays no role in our subsequent 
arguments. 

In many cases of physical interest, the system can be described by averaged diffusion 
coefficients i3)'i'2 - «n -- Once more, the averaging process over TM. for bubble coUisions 
justifies this fact. Then, we can replace the different values of the tensor components by a 
family of scalar functions ^ = 1, 2, . . .} in each term of (9). A class of examples of 

vacua very frequently encountered in field theoretical applications for which such a reduc- 
tion is possible is provided by homogeneous spaces [19] of compact Lie groups. Moreover, 
we make the additional assumption that the function 1^(0), to which the coefficient of (8) 
reduces, varies slowly over M.. This allows us to ignore the contribution of the derivatives 
of D and (8) reduces to the diffusion equation 



where denotes the Laplace-Beltrami operator associated with the metric of A4. Ex- 
amples of a slow variation of D[(f)) are A4 which are symmetric spaces [19]. Even more 
specifically, when is a compact quotient of a sphere or a torus by the free action of 
a discrete group, as in the case of projective spaces for instance, then D{(f)) is exactly 
constant [22]. 

The solution of (11) with the initial condition 



is the heat kernel of the Laplace-Beltrami operator on Ai. Using the maximum 
principle one can prove that P{(p,t), which is really the heat kernel of on Ai, is 



dt 



(11) 



P(0,O) = 5(0) 



(12) 
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positive definite [21] , [23] . Then its normalized form can be legitimately called a probability 
density function, as we have been assuming all along. We are interested in the behavior 
of the system for relatively short times t. Indeed, at very long t the effects of mixing will 
have become significant and, as was remarked after (1), this would introduce additional 
complications that could mask or even invalidate this line of argument. At small t, i.e. 
as t — > 0"*" the heat kernel has the asymptotic expansion [21], [23] 

P(0) = (4.t)-^exp (-J^) Vidm i^tu.m + oit'')^ (13) 

where r}{d{(t))) is a bump function which has the value one, if (j) is in the cut-locus of the 
origin O, and zero if (j) is outside it. We can then clearly see [21] that the leading behavior 
of P(0) as i ^ 0+, is 

lim HnP(0) = -M^ (14) 
which is actually true for any e (Varadhan's asymptotic formula) [21]. 

Let the initial and final values of on be fixed. Then, according to (14), the 
stochastic process described by (2) and (3) gives rise to a probability density P{(j),t) 
which is maximized if the distance function d(0) between O and (p is minimized, i.e. 
along the minimal geodesic joining them. This observation is the geodesic rule. 



III. ON THE APPLICABILITY OF THE GEODESIC RULE 

Let Imin and / be the lengths of the minimal geodesic and some other piecewise smooth 
curve of M. joining O and 0, respectively, and let Pmini'P) ^^'^ P{4') the corresponding 
probability densities at in the limit i — > 0"*". Then (13) implies that 

P(0) = Pmm{<t>) {-^^T^} (15) 

so the contribution of non-minimal paths is exponentially suppressed in expressions like 
path integrals. It is worth noticing that the geodesic rule can actually be violated. We 
have argued that it is not really a rule, but rather a minimalist prescription for calculat- 
ing the density of global topological defects. If someone is interested in calculating the 
number (density) of topological defects with some degree of accuracy beyond a simple 
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estimate, either analytically or numerically, the geodesic rule will provide a fast and accu- 
rate way most of the times. However, we expect deviations from its predictions to become 
important in cases in which there are many paths for which P{4)) ~ Pmin- One such case 
occurs, if, for instance Al is a C; manifold [19], i.e. a compact Riemannian manifold 
all the closed geodesies of which have length /. Examples of Q manifolds are symmetric 
spaces of positive sectional curvature, manifolds which occur frequently as vacua of phys- 
ical systems. 

What we have effectively done in this approach, is that we have arrived at an integra- 
tion measure on the appropriate space of functions on A4, namely the Wiener measure 
on this space [24] . This is then used to perform the thermal path integral describing the 
phase transition which gives rise to the topological defects. The path integral is essentially 
determined by the Feynman-Kac formula [24] , [25] 



for the kernel of e^*"^, where Ti. denotes the effective Hamiltonian on M. of the system of 
colliding bubbles. Then, in the normal coordinate system that we are using and disregard- 
ing any curvature terms, which provide second order corrections in t and arc eventually 
factorized in a path integral, V{t) is trivial on the path space of Ai. Therefore, the dom- 
inant contribution in the thermal path integral comes from evaluation of the Laplacian 
on the shortest possible paths on a fact which establishes the geodesic rule. 

The result of this analysis relies on establishing, for this system, the probability dis- 
tribution given by (13). This probability distribution is Gaussian with respect to the 
distance function d{(j)). The Central Limit theorem [26] states that under some very gen- 
eral conditions a probability density on Ai would converge to (13). Wc arc then led to 
believe that although (13) is derived under some rather strong assumptions, it should be 
valid even in cases that violate one or more of these assumptions, as long as the assump- 
tions of the Central Limit theorem arc obeyed. The geodesic rule describes the leading, 
universal, behavior of such systems. One class of systems that violate the geodesic rule 
are ones whose bubble collisions arc described by a non-Markovian process such as, for 
instance, a Levy process [17]. Even in such cases, variations of the Generalized Cen- 
tral Limit theorem still hold, a fact however which docs not guarantee the validity of the 




(16) 
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geodesic rule in its original formulation, which is the one that we have proved in this paper. 

A case in which it is known that the geodesic rule does not hold is when the dynamics of 
the order parameter is dominated by fluctuations [4], [8], [27]. Such an observation however 
does not contradict our conclusions. Indeed, when fluctuations dominate, the absolute 
minimum of the effective potential [10] is what determines M.. The loop expansion of the 
effective potential, in powers of or a corresponding thermal parameter [10], shows that 
fluctuations force the effective potential to develop singularities [10]. The nature of these 
singularities is not clear, but it becomes amply evident through such a procedure that A4 
can no longer be considered to be a Lie group, a homogeneous space or even a manifold 
as in the classical/non-thermal case [22]. As a result, all the arguments leading to (11) 
and subsequently to (13) and (15), which rely on being a Riemannian manifold do not 
hold any more. We conclude then that the treatment of cases of dominant fluctuations 
lies outside the scope of our approach, hence there is no contradiction with our results. 



IV. DISCUSSION AND CONCLUSIONS 

In the case of a first-order phase transition [11] there is an energy barrier that has 
to overcome in order to move from the meta-stable to the stable vacuum. If we consider 
the Fourier transform of the Fokker-Planck equation, then we expect that D{k) will have 
a strong peak on wave-numbers k comparable to the width of the potential well around 
the meta-stable vacuum. In the case of continuous transitions (spinodal decomposition), 
there is no longer an energy barrier that would act as a localizing factor and would force 
D{k) to have sharp peaks. The arguments presented above do not use this localization 
in momentum space, therefore they are equally valid in both cases. This, of course, does 
not mean that the two cases are physically identical. In a continuous transition, we can 
still consider a region that shares the same value of the scalar field, but does not possess 
a physical boundary, like a bubble wall, as in the case of first order phase transitions. 

To conclude, we have shown that the geodesic rule holds for a general class of systems 
undergoing a phase transition, as a leading approximation. Violations of the rule can 
and will generally occur. The validity of this rule depends on establishing a Gaussian 
process through which the phase transition from the meta-stable to the stable vacuum 
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proceeds. For any system that can be described by such a stochastic process, which obeys 
the assumptions of the Central Limit theorem, the density of the resulting global defects 
can be accurately estimated by using the geodesic rule. 

Acknowledgements: We are grateful to Dr. P. Benetatos for his constructive criticism 
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for several comments that helped improve the presentation of this paper. 
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